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摘 要 : 基于 CFD 方法 ,采用 混合 两 相模 型 、 标 准 kK - s 方程 和 sighal 发 展 的 全 空 化 模型 ， 对 具有 特殊 结构 形式 的 侧 壁 式 压 水 
室 离心 泵 空 化 特征 进行 数值 计算 。 预 测 了 泵 在 未 发 生 空 化 和 空 化 条 件 下 的 能 量 性 能 ， 并 且 分 析 了 设计 流量 下 叶片 表面 和 流 道内 
的 空 泡 分 布 情况 。 结 果 表 明 : 在 大 流量 工 况 下 ， 泵 的 必须 汽 刨 余 量 大 ， 容 易 产生 空 化 现象 。 随 着 泵 进口 压力 的 降低 ， 空 泡 最 先 
出 现在 泵 叶片 背面 靠近 进口 边 附近 ; 当 进 口 压力 进一步 降低 时 ， 空 泡 区 域 开始 快速 增加 ， 向 叶片 出 口 延伸 ， 泵 的 扬程 开始 大 幅 
下 降 。 进 口 压力 继续 降低 时 ， 叶 片 工 作 面 开始 出 现 空 泡 ， 叶 轮 周 疝 流 道内 开始 充满 空 泡 ， 此 时 泵 的 扬程 急剧 下 降 ， 泵 的 稳定 运 
行 条 件 被 破坏 。 不 同 汽 蚀 余 量 下 ， 叶 频 处 压力 幅 值 变化 明显 ， 在 临界 汽 蚀 点 处 几乎 达到 极 小 值 。 
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Numerical Simulation on Cavitation Flow in a Centrifugal 


Pump with Slope Volute 
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Abstract: Full cavitation model、 mixture and standardx - s model were used to investigate cavitation flow in a 
centrifugal pump with slope volute. Performance of model pump at cavitation and non-cavitation conditions were 
obtained, and cavitation bubbles distribution in blade surface and channel were discussed. The results show that: at 
high flow rate, cavitation easily occurs in model pump. When static pressure at pump inlet decreases, cavitation 
bubbles firstly appear in suction pressure side of blade inlet. With pressure reducing further, cavitation region 
expands to impeller outlet, and head of model pump is affected obviously. Then cavitation regions transforms to 
blade pressure side, and the whole impeller channel are filled with cavitation bubbles leading to pump performance 
dropping rapidly. At full cavitation condition, steady operating status of model pump is destroyed. Under different 
NPSH working conditions, pressure amplitude at blade passing frequency differs, and pressure amplitude almost 
reaches minimum value at critical NPSH. 
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剧 下 降 ， 并 且 伴 随 大 幅 振 动 噪声 的 出 现 。 因 此 如 何 


0 二 号 最 大 限度 的 抑制 空 化 现象 的 产生 并 且 准 确 预测 空 
水 力 机 械 中 , 由 于 局 部 压力 的 降低 诱发 的 空 化 化 发 生 位 置 成 为 急 待 解决 的 问题 B` 4。 
现象 一 直 是 研究 的 重 难点 问题 , 空 化 产生 的 振动 品 为 降低 离心 泵 振动 水 平 , 提出 一 种 压 水 室 扩散 
声 和 对 叶片 表面 的 破坏 作用 影响 水 力 机 械 的 安全 段 呈 一 定 倾斜 角 的 侧 壁 式 压 水 室 。 空 化 是 离心 泵 振 
稳定 运行 目 。 空 化 现象 属于 典型 的 非 稳 态 气 液 两 相 动 重要 激励 源 ,因此 深入 研究 侧 壁 式 压 水 室 离心 泵 
运动 ， 空 化 产生 、 发 展 、 省 灭 过 程 涉及 传 质 传 热 ， 内 部 空 化 发 展演 变 规 律 对 认 清 空 化 诱发 振动 关系 
属于 多 学 科 交叉 领域 中。 有 较 大 的 意义 。 
离心 泵 广泛 应 用 于 工业 运输 、 城 市 生活 、 农 业 采用 数值 计算 方法 对 比 转速 130 的 侧 壁 式 压 水 
灌溉 等 领域 。 受 空 化 现象 的 影响 ， 离 心 泵 扬程 将 急 室 离心 泵 进行 分 析 。 基 于 fluent 软件 ， 采 用 sighal 
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建立 的 全 空 化 模型 和 两 相 混合 模 型 在 不 同 工 况 下 对 模型 泵 进行 定常 和 非 定常 空 化 数值 计算 , 得 到 了 
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叶片 表面 空 化 产生 部 位 和 发 展 过 程 规律 以 及 泵 的 
能 量 性 能 变化 趋势 。 
1 离心 稍 基 本 参数 


设计 一 台 侧 壁 式 压 水 室 离 心 条 ， 表 1 给 出 了 模 
型 泵 的 基本 性 能 参数 。 


表 1 模 型 泵 性 能 参数 
Tab. 1 Parameters of model pump 
参数 数值 
流量 0 48 m/h 
扬程 万 7.8m 
转速 n 1450 r/min 
比 转速 六 130 
叶轮 进口 直径 亡 100 mm 
叶轮 外 径 DD， 172 mm 


2 数值 计算 方法 

根据 离心 泵 流动 特性 ， 将 流动 区 域 分 为 4 个 区 
域 静止 区 一 进口 段 、 出 口 段 、 压 水 室 ; 旋转 区 一 
叶轮 。 由 于 空 化 现象 的 复杂 性 , 空 化 区 密度 梯度 大 ， 
低 质量 网 格 容 易 造 成 计算 过 程 发 散 。 采用 ICEM 对 
几何 模型 进行 结构 网 格 划分 , 壁面 附近 进行 加 密 处 
理 , 并 且 进 行 网 格 无 关 性 检查 以 保证 计算 结果 的 精 
度 ， 网 格 总 数 为 120 万 。 


图 1 计算 域 


Fig.1 Computational domains 


图 2 单 流 道 叶轮 网 格 图 
Fig.2 Mesh of impeller single channel 

2. 1 控制 方程 

混合 流体 相连 续 方程 
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Prv(pr) =0 (1) 
空 泡 连续 方程 

TPN +V(pfW) -RR C2) 
动量 方程 

FPD+V (pW) = Vp+3 Vt pv] (3) 


+V[(K+H)V:vIt+pg 
标准 x 一 2 方程 


Be (4) 
如 
SU RE | PE (5) 
Ox, Ox Oo. Ox, 
2 
oO(peu) 6 有 6c Ces G. pe (6) 
Ox, Ox, Oo, /Ox, K K 
Gp 9 
Ox, Ox, JOx; 


Cs=1.44 C,,=1.92 C,=0.090,=1.3 ae=1.0 
0 为 气 / 液 混合 流体 质量 密度 ,7 为 混合 流体 速 
度 ，f 为 空 泡 相 的 质量 组 分 ，Ae 为 水 蒸气 凝结 率 ， 
六 为 静 压 力 ， 4 为 分 子 粘 性 系数 ， 4, 为 沙 流 粘性 
系数 。 
计算 中 的 水 的 温度 设 为 30C, 水 的 饱和 蒸汽 压 
PF=3540Pa, 水 / 空 泡 表面 张力 0=0. 0717 NMm。 假 定 
水 中 不 凝结 性 气体 的 质量 组 分 £=0. 8X10”。 计算 
中 流 场 的 压力 均 为 绝对 压力 。 


2. 2 计算 方法 及 边界 条 件 


固体 壁面 采用 无 滑 移 条 件 , 采用 标准 壁面 函数 
处 理 近 壁 区 域 复杂 流动 。 采 用 一 阶 迎风 格式 对 速 
度 、 消 动能 耗 散 率 及 滑动 能 进行 离散 。 采 用 SIMPLE 
算法 实现 速度 和 压力 之 间 的 耦合 。 在 计算 域 的 进口 
断面 采用 绝对 速度 进口 ， 并 且 在 进口 给 定 x 一 & 的 
数值 ， 该 值 由 式 8) 确定。 采用 压力 出 口 条 件 ， 
并 降低 松弛 因子 以 保证 计算 过 程 稳定 。 


和 


k=0.005u,,,e =(C, rx?) /1 (8) 
1=0.07D, 
3 计算 结果 分 析 
3. 1y+ 检 验 

标准 x 一 & 模型 适用 于 高 雷诺 数 消 流 流动 ， 对 
于 近 壁 面 附近 复杂 层 流 及 过 渡 流 动 状态 ， 标 准 
Kk 一 & 则 无 法 求解 ， 近 壁 区 的 流动 则 采用 壁面 函数 


来 处 理 , 因此 近 壁 面 第 一 层 网 格 节点 应 处 于 满 流 边 [ 加 
界 层 内 。y+ 取 值 范 围 一 般 为 30 一 200。 图 3 给 出 了 sn te Jos 
叶片 表面 的 y+ 值 分 布 ,可 以 看 出 其 变化 范围 满足 要 ,| a 
求 ， 网 格 满足 计算 精度 要 求 。 | / pp ]24] 1ouo 
160 上 6 上 
| | of - 、 i 
120 上 上 学 
+ 100 上 Ms dE He Me 
~ ao[ i 
| 图 4 泵 性 能 曲线 
人 Fig.4 Performance of model pump 
| 下 泵 的 必须 汽 蚀 余 量 NPSHc 分 别 为 1.12 m、1.25 
“| 0 星 时 m、1.7 m。 从 NPSHe 可 以 得 到 大 流量 下 泵 抗 汽 蚀 
CC O10 lB 性 较 差 ， 这 与 理论 分 析 是 一 致 的 。 随 着 进口 压力 的 
到 3 叶片 表面 y+ 分 布 进一步 降低 ， 泵 的 扬程 会 急剧 下 降 ， 严 重 影响 离心 
Fig.3 y+ distribution of impeller 泵 的 稳定 运行 。 
3. 2 侧 壁 式 压 水 室 离心 泵 能 量 性 能 分 析 sl 0 
人 
图 4 给 出 了 侧 壁 式 压 水 室 离心 硝 性 能 曲线 , 友 ER 
的 最 高 效率 点 偏 到 1. 05By 附 近 , 在 设计 流量 点 上 ， Ei 1 
泵 的 扬程 为 了 4m, 基本 满足 设计 要 求 。 | 有 
图 5 给 出 了 侧 壁 式 压 水 室 离心 泵 3 个 典型 工 况 | 六 
下 (0. 8@By, 1.0Cw 1. 2Bw) 的 空 化 性 能 曲线 。 当 进 
口 总 压 较 大 时 ， 泵 内 无 空 化 现象 产生 ， 离 心 和 泵 能 量 LE RE 
性 能 未 受 影 响 ， 扬 程 几乎 保持 不 变 。 随 着 进口 总 压 2 


5 离心 泵 空 化 性 能 曲线 
Fig.5 Cavitation performance of pump 


的 降低 ， 离 心 泵 叶片 表面 开始 产生 空 泡 ， 使 得 泵 扬 
程 开始 下 降 。 根 据 扬程 下 降 3% 准 则 可 得 3 种 工 况 


体积 分 数 “ 国 二 二 二 


RAR A 


NPSH=5. 15 mm NPSH=3. 18 m NPSH=2. 66 m NPSH=2. 19 m 
NPSH=1. 69 m NPSH=1. 25 m NPSH=0. 96 m NPSH=0. 92 m 
图 6 叶片 背面 空 泡 分 布 


Fig.6 Bubbles distribution on blade suction side 
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片 表面 空 泡 体积 分 布 图 


。 从 图 中 可 以 得 


化 时 ,叶片 进口 边 是 低 
时 ， 低 于 液体 汽化 压力 ， 


压 区 ， 因 此 当 进 


图 6 给 出 了 1.0Cx 工 况 下 ， 不 同 NPSH 值 时 叶 


到 : 在 无 空 
总 压 降低 


空 泡 将 最 先 在 


叶片 进口 边 


附近 产生 ,此 时 空 泡 区 域 和 空 泡 体 积分 数 还 处 于 较 
低 水 平 ， 空 化 对 泵 的 性 能 影响 较 小 。 当 压力 进一步 
降低 时 ,最 大 空 泡 体 积分 数 区 域 还 处 于 叶片 进口 边 
附近 ， 空 泡 区 域 开始 增加 ， 向 叶片 出 口 边 延 促 。 此 
时 泵 扬程 受到 影响 ， 当 NPSH=1.25 m 时 ， 泵 扬程 
下 降 3%。 进 一 步 降低 进口 压力 时 ， 空 化 区 域 继续 
增加 ， 空 泡 运动 到 高 压 区 时 产生 涡 灭 ， 形 成 明显 的 
气 / 液 分 界面 。 当 NPSH=0.92 m 时 ， 空 化 区 域 占据 
整个 叶片 1/2 左右 ， 此 时 泵 的 扬程 急剧 下 降 ， 工 作 
性 能 恶化 , 空 泡 体 积分 数 最 大 值 出 现在 叶片 前 盖 板 
附近 


3 空 化 流动 非 定常 计算 结果 


离心 汞 叶轮 内 部 产生 空 化 现象 时 , 稳定 的 流动 
结构 被 破坏 ， 空 泡 初生 、 发 展 和 溃 灭 过 程 将 对 叶轮 
内 流动 产生 明显 影响 。 叶轮 一 隔 舌 动静 干涉 作用 诱 
发 的 压力 脉动 是 诱导 离心 泵 低频 振动 的 重要 激励 
源 ,。 为 了 得 到 空 化 流动 对 泵 内 部 压力 脉动 的 影响 规 
律 , 在 1.0Bw 工 况 下 对 侧 壁 式 压 水 室 离心 泵 进行 空 
化 非 定 常数 值 计算 ， 时 间 步 长 设置 为 1.15x10 s。 
监测 点 设置 如 图 7 所 示 。 


焉 


图 7 监测 点 
Fig.7 Monitoring points 


图 8 给 出 了 3 个 不 同 监 测 点 处 , 叶 频 处 压力 脉 


动 幅 值 随 有 效 空 化 余 量 变化 曲线 。 开 始 阶段 ， 叶 片 
进口 压力 远 高 于 汽化 压力 ， 无 空 化 现象 产生 ， 压 力 
脉动 幅 值 较 大 。 随 着 叶轮 进口 压力 的 降低 ， 有 效 空 
化 余 量 减 小 ,压力 脉动 幅 值 呈 下 降 趋 势 ， 此 时 模型 
的 性 能 几乎 不 受 影响 。 当 达到 临界 汽 蚀 余 量 点 时 ， 
叶片 进口 边 空 泡 区 域 扩 大 , 此 时 压力 脉动 基本 达到 
最 小 值 。 有 效 空 化 余 量 进一步 降低 时 ， 空 泡 充 满 叶 
轮流 道 ， 此 时 泵 的 流动 状态 完全 被 破坏 ， 导 臻 压力 
脉动 幅 值 迅速 上 升 。 因 此 可 以 将 压力 脉动 迅速 上 升 
的 拐点 作为 判断 泵 内 临界 汽 蚀 状态 的 依据 。 
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四 


250 一 as 一 p0 一 * 一 p2 一 人- p6 


E 
< 
NPSH (m) 
图 8 叶 频 幅 值 变化 规律 
Fig.8 Pressure pulsation amplitude 
4 结论 


应 用 数值 计算 方法 对 侧 壁 式 压 水 室 离心 泵 进 
行 空 化 数值 分 析 ， 得 到 以 下 结论 : 
1、 大 流量 下 ， 泵 的 必需 汽 刨 余 量 大 ， 泵 容易 产生 


因 


Ds 空 泡 最 先 出 现在 叶片 年 进口 边 附近 ， 进 口 压 
压 区 溃 灭 形成 明显 的 汽 \ 液 分 界面 ， 压 力 继续 降低 
时 ， 流 道内 充满 空 泡 ， 泵 的 扬程 急剧 降低 。 

3、 临 界 汽 蚀 点 处 ， 叶 频 处 压力 脉动 幅 值 几乎 达到 
最 小 值 , 之 后 随 着 有 效 汽 蚀 余 量 的 降低 压力 幅 值 迅 
速 上 升 。 
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